Spinon signatures in the critical phase of the (1, ~) ferrimagnet in a magnetic field 
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We propose an effective theory for the critical phase of a quantum ferrimagnetic chain with alter- 
nating spins 1 and | in an external magnetic field. With the help of the matrix product variational 
approach, the system is mapped to a spin-i XXZ chain in an (effective) magnetic field; as a byprod- 
uct, we obtain an excellent description of the optical magnon branch in the gapped phase. Recent 
finite-temperature DMRG results for the low-temperature part of the specific heat are well described 
by the present approach, and the "pop-up" peaks, developing near the critical field values and in the 
middle of the critical phase, are identified with the contributions from two different spinon bands of 
the effective spin-| chain. The effect should be as well observable in other spin-gap systems in an 
external field, particularly in spin ladders. 



75.10.Jm, 75.50.G, 75.40.Cx, 75.40.Mg 

Introduction. Recently, there has been a considerable 
interest in the properties of one-dimensional (ID) quan- 
tum fcrrimagnets.LTa A generic example of a ID ferrimag- 
net is the Heisenberg spin chain with antiferromagnetic 
nearest neighbor interaction and alternating spins 1 and 
5, described by the Hamiltonian 

H = jJ2(S n r n + T n S n+1 ) -h^iSl+O, (1) 

n n 

where S n and r„ are respectively spin-1 and spin-i op- 
erators in the n-th elementary magnetic cell, and h — 
g^sH j where H is the external magnetic field. According 
to the Lieb-Mattis theorem,E2l for h = the ground state 
of the system has total spin S'tot — L/2, where L is the 
number of unit cells, and thus is necessarily long-range 
ordered, making the probltgu. amenable to the spin wave 
theory (SWT) approach .Mia Since the elementary cell 
consists of two spins, SWT yields two types of magnons: 
a gapless "acoustical," or "ferromagnetic" branch with 
S z = L/2 — 1, and a gapped "optical," or "antiferromag- 
netic" branch with S z = L42 + 1. The optical magnon 
gap was numerically founds to be A opt = 1.759 J. The 
existence of two magnon brancbes,manifests itself in var- 
ious thermodynamic quantities H : Q 

In magnetic field the acoustic branch acquires a gap, 
while the optical gap decreases with the field. If the field 
exceeds the critical value h = h c i = A opt , the optical gap 
closes, and the system enters the critical phase, which is 
expected to be of the Luttinger-liquicL-tjype, in analogy 
with other models jl<L£. , spin-1 chainjiil gapped spin-i 
chains and laddersEj'Ej) . This phase extends up to the 
second critical field h = h C 2 = 3 J, where there is an=. 
other transition to the saturated ferromagnetic phase.cl 



In the gapless regime h c i < h < h c2 , the tempera- 
ture dependence of the specific heat has revealed a puz- 
zling behaviorS with a single well-pronounced low-T peak 
which pops up around T rs 0.2 J when h is in the middle 
between h c \ and h C 2] when h is shifted towards h c \ or 
h c 2, the peak becomes flat and develops a shoulder with 
another ill-pronounced peak centered at about 0.05 J. 

In this paper, we show that in the critical phase the 
model (|l|) can be mapped to an effective spin-i XXZ 
chain in an external field, which yields quite an accu- 
rate description of all the features of the low-temperature 
part of the specific heat. The "pop- up" peak structure 
described above can be explained by the presence of two 
different spinon bands of the effective spin-i chain. This 
effect is rather general and should be observable in other 
gapped ID spin models, e.g., spin-i ladders. 

Effective model. The general idea of any mapping 
to an effective model is to reduce somehow the Hilbert 
space of the problem, keeping only a few "most impor- 
tant" states per the elementary cell. For example, in the 
strong-coupling limit of the spin-i ladder one keeps for 
each rung only the singlet and the lowest-energy triplet 
component J13 For the (1, 5) ferrimagnet the set of the 
cell wave functions ipj m consists of a doublet {j = \) and 
a quartet (j — |). One would naively expect that in a 
strong field it is now necessary to keep three states with 
(j,m) = (|,|), (§, |), (§, |); it is clear, however, that 
there is only one ground state and one excitation becom- 
ing gapless at h = h c i, so that the interplay should be 
effectively between two cell wavefunctions being "proper" 
linear combinations of ipj m . Thus, the problem is how to 
identify those proper combinations. The first step is to 
describe accurately the optical magnons, since the criti- 
cal phase is formed by their condensation into the ground 



1 



state. 

The linear SWT captures tbe essential physics of the 
model at the qualitative level;cl however, for a quantita- 
tive description one has to go far beyond the linear 
proximation, keeping the higher-order terms in 1/S 
e.g., in linear SWT A opt =i-j/, and the corrections from 
the leading 1/5 term yielcfl 1.676 J instead of the nu- 
merical value 1.759 J. Thus we use a different scheme 
which has proveofi to be very successful for the (lwil 
ferrimagnet, namely the variational matrix productO't£l 
(MP) approach. In Ref. ||the following (non-normalized) 
variational ground state wave function was proposed: 

*o = tr( 5 i.g 2 ---. 9L ), (2) 

M 1 '* = — l —a°ibx i + \ -<r +1 ipi i, 
x/3 2 ' 2 V3 ^3-2' 

3 1 1 1 

M 1( 5 = —=o~ x %bz s H =cr +1 ^3 i =a°ib3 i . 

V2 2 ' 2 V6 21 2 V3 2,2 



-ifTy) are the Pauli ma- 



Here fj° = er z and (T ±:L = ^-j=(a x 

trices, 1 is the 2x2 unit matrix, ipj m denote spin states 
of the n-th cell, and it, v are the variational parame- 
ters whose optimal values minimizing the energy of (|l|) 
at h = were foundB to be u = 1.3026, v = 1.0788, 
with the variational ground state energy per unit cell 
E™ — —1.449 J (to compare with the numerical value 
E n = —1.455). It was shown that \&o possesses the cor- 
rect quantum numbers of S^ot = S* ot = L/2 and provides 
a very good description of the ground state correlations. 
Any quantum averages for MP wave functions are readily 
calculated using the transfer matrix technique; this nice 
formalism was first proposed in Rcfs. [TJJl^ The MP ap- 
proach is especially well suited to this problem since the 
fluctuations are extremely short-ranged, with tJae corre- 
lation radius smaller than one unit cell length.ciu 

Now we construct the simplest trial wave function for 
the optical magnon with the momentum k: 



\k) =2^e* n \n), \n) = tr (g x ■ ■ ■ 9n-i9n9n+x ' • -Ql), 

n 

9n = 1^4,3 +wcr +1 Vi,i + VWS/Q 1 ' 2 , (3) 
The form of g n in (J3j) is dictated by the requirement that 



is the state with St, 



cording to the general formalism presented in Ref. a g n 
should carry the "hyperspin" quantum numbers (3, §), 
while each g n carries (A, | ). Here /, w are still free (vari- 
ational) parameters. Note that generally the states |n) 
are not orthogonal to each other, but are orthogonal to 
^o- For the following, it is important to note that g n can 
be represented in the form 



9n 



f - 1 +1 



V2 



<? +1 9n 



w er +1 tpi i 



so that |fe) in fact depends only oihd= {w+\/2(u+ . 
Thus, one parameter in (|3|) is redundant, and we fix that 
by choosing 



f={z + wV6(v - uV3)}/(2 + z) , 
z ee 2\/3 uv + {3 + 12v 2 (u 2 + 1)} 1/2 . 



(4) 



which yields a remarkable property: the set of one- 
magnon states {\n}} becomes mutually orthogonal, con- 
siderably simplifying further calculations. Moreover, one 
can show that arbitrary 7V-magnon states \n\, . . . njv) be- 
come orthogonal; however, they are not normalized, and, 
generally, (n, n + m\n, n + m) does not coincide with 
(n\n) 2 , though it tends very fast to the latter value as 
m increases. The one-magnon norm and the matrix ele- 
ments of the Hamiltonian (Q) are given by 

Wo ee (n|n) = 0.1173w 2 - 0.0337k; + 0.1432 , 
(n\H\n') = J X \ n - n '\A + BS ]n _ n , Ul + C5 nn> ) , 
x = 0.2147, A = -0.0483w 2 + 0.0808u> - 0.0284 , 
B = -0.0839W 2 + 0.0379k; + 0.0010 , (5) 
C = 0.3359w 2 - 0.0495w + 0.5860 , 



where H = H — LE™ 1 ' . Now one can calculate the 
dispersion by minimizing the excitation energy e(k) = 
{k\H\k)l{k\k}, 



C + 2Bx cos k + A 



l-x 2 



1 + x 2 — 2x cos k 



}• (6) 



In principle, the optimal value of w should be calculated 
separately for each k. One might first try to minimize the 
gap A opt = e{k = 0), which yields w — w — —3.8605, 
with A™t ~ 1.754 J, the resulting dispersion for w(k) = 
wo being in excellent agreement with the exact diagonal- 
ization data (see Fig. 

Now, we can obtain the desired mapping by intro- 
ducing effective spin-i states \a n ) at each cell, a n = 
±5, and making the identification \a\a.2 • • • oil) — 
tt(RiR 2 ■ ■ ■ Rl), where R n = g n , g n for a n = \, -3, 
respectively. Introducing the spin-^ operators s^' , and 
restricting all effective interactions to nearest neighbors 
only, one can write down the effective Hamiltonian for 
the critical phase of the (1, 4) ferrimagnet (111): 



L/2 + 1, then, ac- jj e _ 



_ \ " Jxy , + 



^ 2 



( s n s n+l + s n s n+l) + ^z S n S n+l 2/l e S n , 



J xy = -2t u J Z = U U h B = (h-t -U 1 )/2, (7) 

_ (n\H\n + m) _ (n,n + m\H\n,n + m) 
I'm. — ~ , yJ 7n — ~i ' 1 ' r 2£q . 



(n, n + m\n, n + m) 



Here a remark is in order. From (|5|) one can see that the 
hopping amplitudes t m cx x m are very small for m > 2, 
t 2 = -0.0242, t 3 = -0.0052, . . ., which justifies keeping 
the terms of the s+s~ +m type only up to m = 1 in (Q). 
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The same reasoning holds for neglecting other interaction 

terms like (s+ s n+l s n-l s n+2)> (s+ s n+2 s n s n+2): etc - 

The numerical values for the parameters to, ti, U\, 
which determine our effective model, are to = 2.3370, 
ti = -0.2608, Ui = 0.1187, so that the spin-i chain (g) 
is in a gapless XX phase with A = J z /J xy ~ 0.227. At 
the critical valueEO of the field h e = ft 6iC = hJ xy {l + A) 
the spin-i chain undergoes a transition into the saturated 
ferromagnetic phase. The obvious symmetry h e t— > — h e 
corresponds for the ferrimagnet to the symmetry around 
ho = \ih c i + h C 2). Such a symmetry was indeed observed 
in Ref. Q for the behavior of the specific heat, though for 
the h < ho side it was somewhat plagued by the low-T 
contribution of the acoustical magnon branch. In what 
follows, we compare our results only with the data for 
h > ho. Because of cutting out non- nearest neighbor in- 
teractions, the critical fields for the model ([?]) are slightly 
different from the real ones: h e = ±/i e c corresponds to 
ft = h' cl ~ 1.81 and h = h' c2 ~ 3.09 instead of the cor- 
rect values listed above. Thus, for a comparison to the 
numerical data below, we have linearly rescaled the field 
interval {h' cll h' c2 ) onto (h c i,h C 2)- 

The specific heat. The temperature dependence of the 
specific heat for the. spin-i chain can be calculated using 
Kliimper's versionllj of the thermodynamic Bethe ansatz; 
this amounts to solving numerically the following system 
of two coupled nonlinear integral equations: 



yi 



K(3 



ir(3h 



7 cosh 22 2(ir — 7) 
K(x — x' — i"f + i<5)-L(?/2)| 



+00 

dx'\K(x- x')L( yi ) 



(8) 



the other equation being of the same form with i — > — i, 
yi 2/2, and h — > —ft. Here L(y) = ln(l + e~ y ), 
A = cos 7, 6 is an infinitesimal positive number, h — 
4ft e /( J xy sin 7), — 2T/(J xy smj), and the kernel K(x) 
is defined as follows: 



K{x) 



I f 30 ^ sinh[(7r/2- 7 ) g ] cosgx 



-in 



cosh(7g/2) sinh[(7r — j)q/2] 



(Note that the last term in (||) contains a singularity!). 
Then the specific heat C(T) can be found as 



C = f3 



2^1 

9/3 2 ' 



/ = 



20 



+°° L( yi ) + L(y 2 ) 
ax ■ 



cosh(7rx/7) 



(10) 



In Fig. □ we compare the results given by (fl(J with the 
numerical data for the low-temperature part of C (T) for 
the ferrimagnet (Q), obtained by means of the recently 
proposedlia transfer-matrix density matrix renormaliza- 
tion group (DMRG) technique. The agreement is gen- 
erally satisfactory, and becomes very good near h = fto- 
The "pop-up" peak behavior is clearly reproduced. 

Qualitatively, the behavior of C(T) can be understood 
without appealing to the machinery of the Bethe ansatz. 



Passing to the Jordan- Wigner fermions in (j7j) , and apply- 
ing the first-order perturbation theory in A, one obtains 
the following simple expression for the fermion dispersion 
e(k): 



e(k) = A — A — cos k 

- (2A/tt)0(1 - A){arccos 



(11) 

A + (l-A 2 ) 1/2 cos/fc}, 



here e{k) is measured in J xy units, A = 2h e /J xy , and 9 
is the Heaviside function. The specific heat can be triv- 
ially calculated; in Fig. ||a, we have included two curves 
obtained within the above approximation. One can see 
that, despite the general tendency to overestimate C, 
this simple approach nevertheless captures the essential 
physics. In zero field the contributions into the specific 
heat from particles and holes are equal; with increasing 
field, the hole bandwidth grows up, while the particle 
bandwidth decreases, and the average band energies do 
not coincide. This leads to the presence of two peaks 
in C(T): holes yield a strong, round peak moving to- 
wards higher temperatures with increasing the field, and 
the other peak (due to the particles) is weak, sharp, and 
moves to zero when h e tends to h e ^ c . In the spinon lan- 
guage, one has essentially the same picture for the two 
bands of spipens with opposite spin (the fact first noticed 
by KlumpeniS). 

Summary. An effective theory for the critical phase of 
a quantum ferrimagnet with alternating spins 1 and \ in 
external magnetic field is proposed. Using the matrix- 
product formalism, we have mapped the low-energy sec- 
tor of the system to the spin-i XXZ model in an effec- 
tive magnetic field; as a byproduct, we obtain an ex- 
cellent description of the optical magnon branch in tha 
gapped phase. Recent transfer matrix DMRG resultstl 
for the specific heat C(T) of the ferrimagnet find natu- 
ral explanation within the present approach. The origin 
of "pop-up" peaks observed near the critical field val- 
ues and in the middle of the critical phase, is clarified: 
they can be identified with the contributions from two 
different spinon bands of the effective spin-i chain. Wc 
believe those results should be accessible to experimen- 
tal verification, and we would like to emphasize that this 
effect is general and should manifest itself in the criti- 
cal regime of essentially any gapped ID spin system in 
a magnetic field, as far as it is possible to construct the 
mapping to a spin-i chain. Spin-i ladder compounds like 
Cu2(CsHi2N2)2Cl4 would be the natural candidates l3E!l 
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FIG. 1. Dispersion of the optical magnon branch (in J 
units). The solid line is variational result (^), and symbols 
denote the Lanczos data. 
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FIG. 2. Temperature dependence of the specific heat for 
the model (|l). Symbols denote the DMRG data, and thick 
curves corres )ond to the Bethe ansatz result ( |lo|) for the effec- 
tive model (|7). Two thin curves in (a) illustrate the first-order 
perturbative result based on (|ll|). 
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